Synthetic magnetic fluxes on the honeycomb lattice 
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We devise experimental schemes able to mimic uniform and staggered magnetic fluxes acting 
on ultracold two-electron atoms, such as ytterbium atoms, propagating in a honeycomb lattice. 
The atoms are first trapped into two independent state-selective triangular lattices and are further 
exposed to a suitable configuration of resonant Raman laser beams. These beams induce hops 
between the two triangular lattices and malce atoms move in a honeycomb lattice. Atoms traveling 
around each unit cell of this honeycomb lattice picic up a nonzero phase. In the uniform case, the 
artificial magnetic flux sustained by each cell can reach about two flux quanta, thereby realizing a 
cold atom analogue of the Harper model with its notorious Hofstadter's butterfly structure. Different 
condensed-matter phenomena such as the relativistic integer and fractional quantum Hall effects, as 
observed in graphene samples, could be targeted with this scheme. 



I. INTRODUCTION 

The abihty to produce and trap dilute degenerate Bose 
and Fermi gases [l|, ushered ultracold atoms as pow- 
erful players to explore phenomena mostly studied until 
now in condensed-matter physics, see 0, for compre- 
hensive reviews. Paradigmatic examples are the obser- 
vation of the Abrikosov vortex lattice [1, Q , the Mott- 
superfluid transition [7|, the BEC-BCS crossover d, H] 
and the Kosterlitz-Thouless transition (lo| . 

One current motive in the field is to address the physics 
of two-dimensional (2D) electrons exposed to a magnetic 
field and the associated integer (IQHE) and fractional 
(FQHE) quantum Hall effects Even though atoms 
are neutral, people have rapidly realized that, because 
of the similarity between the Lorentz and the Coriolis 
forces, repulsively interacting atoms loaded in rotating 
traps could target the fractional Hall effect regime fl^. 
Unfortunately the first attempts failed as it is difficult 
to reach the regime of one vortex per particle (Tsl . [l^ . 
One key reason is that the observation of the Laughlin 
states requires a very fine tuning of the rotation frequency 
compared to the trap frequency. 

In the meantime other promising methods using light- 
induced gauge fields (l5l - [ig | where theoretically explored, 
see [20] for a review. For atoms trapped in 2D optical lat- 
tices, these schemes amount to control the phases of the 
hopping amplitudes such that non-intersecting loops ac- 
quire a nonzero phase. This simply mimic the effect of 
a constant magnetic field leading to Harper's model |t21j 
and Hofstadter's butterfiy More elaborate schemes 
targeting non-Abelian gauge fields have been proposed 
[23., ^4]. From an experimental point of view, the success- 
ful generation of an Abelian gauge field in the bulk has 
been reported with the observation of vortices (25l. [2^. 

The aim of our work is to devise workable experimen- 
tal schemes that would create artificial Abelian gauge 
fields acting on a honeycomb lattice filled with earth- 



alkaline atoms. These artificial gauge fields mimic effec- 
tive magnetic fields imposing uniform or staggered fiuxes 
through the lattice. Related schemes have been theoret- 
ically proposed in [l9l [27| for the rectangular lattice. 
Recently a scheme mimicking an effective periodic mag- 
netic field has been proposed in [28j. Our primary in- 
terest in the honeycomb lattice is that its ground and 
first-excited bands give rise to two Dirac points around 
which the dispersion is linear. When loaded with ultra- 
cold fermions around half-filling, the Fermi energy slices 
the band structure around the Dirac points and one gets 
a cold atom analogue of graphene [2^ [s^l- Key exper- 
iments with graphene samples have revealed a particu- 
lar quantum Hall effect [sil Is^ . This is due to the rel- 
ativistic nature of the low-energy electronic excitations 
which behave like massless Dirac particles [s^. Recently 
the first observations of a FQHE state a.t v = ±1/3 in 
suspended graphene were reported ||33, Because of 
the richness and flexibility of the cold atoms technology, 
we believe that experiments where fermions, bosons or 
fermion-boson mixtures are loaded in the honeycomb op- 
tical lattice and are exposed to an artificial magnetic field 
yield situations difficult to explore in graphene research, 
in particular when the hopping amplitudes are imbal- 
anced [s^. It is worth mentioning that configurations 
giving rise to Dirac points in the square lattice have been 
proposed (37l.l38j. However it can be shown that the Dirac 
points are robust in the honeycomb optical lattice where 
the laser beams need not be at perfect respective in-plane 
120° angles from each other or carrying exactly the same 
intensities [sO]. In this context, it is worth noticing that 
the first experimental realization of spin-state dependent 
lattices realizing a honeycomb lattice and loaded with 
ultracold ®^Rb ground-state atoms has been recently re- 
ported [39]. 

The paper is organized as follows. In Section |ll] we 
introduce the laser configuration producing the two un- 
coupled state-dependent triangular lattices we need. In 



2 



Section Iml we introduce the Raman laser configuration 
which induce hops between the two previous sublattices, 
thereby reahzing a (shghtly non regular) honeycomb lat- 
tice where the atoms can move. We show that a single 
Raman laser scheme is not sufficient to induce an artifi- 
cial magnetic field and we propose a 4-beam scheme to 
achieve a uniform flux through the lattice. About two flux 
quanta per unit cell can be generated with our scheme, 
thereby reaching the strong field limit and realizing the 
cold atom analogue of the Harper model. In Section ITVl 
we introduce two set-ups that allow to achieve a stag- 
gered flux through the lattice similar to the one shown 
in |,19i] . About one flux quanta per unit cell can be gen- 
erated with both schemes. In Section |Vl we give the as- 
sociated Harper models. We summarize and conclude in 
Section ED 



II. STATE-DEPENDENT TRIANGULAR 
OPTICAL LATTICES 

In this section we introduce the laser conflguration 
that produces two triangular state-dependent optical lat- 
tices. Like in fl9'|, all our calculations are done for the 
bosonic isotope -'^^'^Yb of ytterbium (iol. I4H. We restrict 
our analysis to two internal states, namely the ground 
state ^5*0, hereafter denoted by \g), and the long-lived 
(lifetime about 20 s) metastable excited state '^Pq, here- 
after denoted by |e). The energy separation between these 
two states is htoo = hcko = hc/Xa, where Ao = 0.578/xm 
and the so-called "magic" and "anti- magic" wavelengths 
of this isotope are Am = 0.76/im and Xam = 1.12/im. 

In a nutshell, our strategy is the following. A first laser 
configuration working at the "magic" frequency creates 
a honeycomb potential with lattice constant a oc Xm 
which confines both Yb internal states in its minima. 
At this stage there is thus no special spatial organiza- 
tion of the Yb internal states among the two triangular 
sublattices of the honeycomb lattice. A second laser con- 
figuration working at the "anti-magic" frequency is then 
superimposed to the previous one to create an optical 
standing- wave potential. When 3a = Xam, as a net re- 
sult of the combination of these two optical potentials, 
state- 1 (7) atoms are trapped in the minima of a triangular 
lattice while state- |e) atoms are trapped in the minima 
of another shifted triangular lattice. As a whole, we get a 
(slightly non regular) honeycomb lattice where state-\g) 
atoms are solely trapped in one of its sublattice while 
the state- 1 e) atoms are trapped in the other. Consider- 
ing the hexagonal Bravais Wigner-Seitz cell of the new 
honeycomb lattice, this means that its vertices are now 
alternatively occupied by state- jg) and state- |e) atoms. 
At this stage however, the state- jg) and state- |e) sublat- 
tices are still uncoupled as atoms trapped in one sublat- 
tice cannot fiip their internal state and hop on the other 
sublattice. How to couple these two sublattices will be 
the topic of Section iHll 

Throughout this Section and Section IIIIl for conve- 



nience purposes, we will use Xamf^ = 0.56/iTO as our 
space unit. 

A. "Magic" honeycomb lattice 

The basic laser configuration creating the honeycomb 
lattice was presented and analyzed in [30]. In the present 
case, we choose a slightly modified version of it by con- 
sidering the superposition of three linearly-polarized run- 
ning monochromatic waves at angular frequency ujm = 
ckm = 27rc/Am with wavevectors — km&i, where: 

§1 = cos6'ey + sin6'ez, (la) 

/I V^^x ^ . • rt ^ /-II \ 

e2 = — cosfc* —+sinOez, (lb) 

83 = cost' —+smOe2, (Ic) 

6 being the elevation angle of the "magic" beams off the 
(Ox, Oy) plane, see Fig. [T] 




FIG. 1. The laser configuration creating a honeycomb optical 
lattice in the {Ox, Oy) plane at the magic wavelength. It is 
made of three monochromatic linearly-polarized beams with 
the same off-plane elevation angle 6. The in-plane projections 
of the laser wavevectors have respective angles 120°. 

Up to an inessential additive constant, the resulting 
optical dipole potential K„(r) is translation- invariant 
along the Oz direction and displays a honeycomb struc- 
ture in the (Ox, Oy) plane with lattice constant a = 
2A„i/(3-\/3 cos6'). After a suitable choice of space and 
time origins, it is then given by V^(r) — Vq[3 + 2/m(r)] 
where: 

/m(r) = cos(bi • r)+cos(b2 • r)+cos ((bi -l-ba) • r)), (2) 

and where Vq denotes the strength of the potential. The 
vectors bi = — ki and b2 = ki — k2 feature the prim- 
itive reciprocal lattice vectors of the honeycomb lattice. 
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In turn they define the primitive honeycomb Bravais lat- 
tice vectors ai and a.2 through • hj = 2TrSij . In terms 
of our space units, we find: 

/m(r) = 2 cos (n^x) cos (V^n^y) + cos {2TT"fx), (3) 

where: 

VSKim cos 9 



2\r. 



(4) 



Oy 




FIG. 2. Hexagonal Wigner-Seitz cell of the honeycomb lattice 
associated to the optical potential generated by Eq. ((3]). The 
figure shows the two primitive Bravais vectors ai, the three 
vectors cj connecting each honeycomb lattice site to its three 
nearest neighbors and the honeycomb lattice constant a. 

Throughout this paper, we assume that a suitable ex- 
ternal confinement along axis Oz restricts the atomic 
dynamics in the (Ox, Oy) plane. Then, for Vb > 0, the 
atoms are trapped in the minima of the potential which 
coincide with its zeroes. These are organized in a regular 
honeycomb structure made of two shifted identical tri- 
angular sublattices [30|. Each lattice site is connected to 
its three nearest neighbors by the vectors (a = 1, 2, 3) 
satisfying: 



ai = ci - C2, 

a2 = Ci - C3, 



(5a) 
(5b) 

(5c) 



Figl^l shows the hexagonal Wigner-Seitz cell of this regu- 
lar honeycomb lattice together with the primitive Bravais 
and nearest-neighbor vectors. The potential is maximum 
at its center, vanishes at each of its vertices and exhibits 
a saddle-point at the mid-point of each of its sides (s^l • 



B. Creating state-dependent triangular sublattices 

The previous "magic" honeycomb potential trap atoms 
in its minima irrespective of their internal state. Start- 
ing from this situation, we would like now to selectively 



trap atoms with a given internal state in a given trian- 
gular sublattice of a new honeycomb lattice. As a result 
the six vertices of the corresponding Wigner-Seitz cell 
would alternate trapped internal states. This is achieved 
by shining two counter-propagating laser beams along 
the Ox direction, their common angular frequency being 
Warn = ckam = 27rc/Aam- Taking the origin of space at a 
point where the "magic" honeycomb potential Kn(r) is 
maximum, this "anti-magic" standing-wave potential is 
given by: 



Krn(r) =^1 



1 -f COS (b • (r - ro)) 



(6) 



where b ~ 2kam&x is the anti-magic reciprocal lattice 
vector, Vi being the corresponding potential strength. 
The position Tq where this potential reaches its maxi- 
mum value is determined by the relative phase of the 
two interfering "anti-magic" laser beams. The effective 
potentials for the ground and excited states then read 
respectively: 



Vgir) = Vrair) + Vamir), 
V;(r) = Kn(r)-K™(r). 



(7a) 
(7b) 



The strategy is now to find parameters for which both 
Vg and Ve sustain the same (regular) triangular Bra- 
vais lattice, both potential minima being organized in 
two shifted triangular lattices. As a whole, one would get 
a honeycomb lattice where adjacent sites would now be 
loaded with atoms having different internal states, or, in 
other words, where each sublattice sustains a given inter- 
nal state. These sublattices of minima will be triangular 
if the "magic" and "anti-magic" lattices match. This the 
case if b belongs to the reciprocal lattice of the "magic" 
honeycomb lattice. This requirement is simply met when 
b = bi + b2, i.e. when Xam — 3a (or equivalently 7 = 1). 
This is implemented by stretching the honeycomb lattice 
and choosing the common polar angle 9 of the "magic" 
laser beams to satisfy: 



cos 6* — 



2An 



V3Aa 



(8) 



For the chosen internal states of the Yb atom, one finds 
the required off-plane angle 9 — 31.41°. 

Next we require a;o — 3a/ 8 (= 1/4 in dimensionless 
units) such that cos(27ra;o) = 0, see Fig. [3l From an ex- 
perimental point of view, this would require to control 
the phases of the laser beams. For this particular choice 
of xq , one has (in dimensionless units) 



Vamir) = Vi 



1 + sin(27ra;) 



(9) 



The state-dependent potentials satisfy Vg{x,y) = Vi + 
Ve{—x,y). Since they are also even functions of y, they 
also satisfy V"g(r) = Vi -I- ye(— r). This means that Vg and 
Ve are thus simply obtained from each other by mere 
reflection about the Oy axis and also by mere inversion 
about the origin. 
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FIG. 3. [Color online] Top: contour plot of the stretched 
"magic" honeycomb potential Vm {a = 2/3 in space units 
of Aam/2). Bottom: plot of the "anti-magic" standing- wave 
potential Vam for a relative position xo — 3a/ 8 = 1/4. The 
"anti-magic" potential takes its mean value half-way between 
the minima of the "magic" potential (vertical dotted lines). 
The superposition of these two potentials gives the state- 
dependent trapping potentials Vg.c = Vm ± Vam- The posi- 
tions of the resulting |(;)-sites and |e)-sites are marked by 
the corresponding blue and red letters. Dash-dotted blue line: 
triangular unit cell of the state-|gi) lattice. Dashed red line: 
triangular unit cell of the state-|e) lattice. The combination 
of these two regular triangular lattices produces a honeycomb 
lattice with alternating state-dependent sites. The unit cell of 
this honeycomb lattice (black solid line) is no longer regular, 
see Fig. \S\ 



Figure H] shows Vg as a function of x when y = 1/ a/3 
for various Vi/Vq ratios. As Vi is increased from zero, the 
two wells of the "magic" honeycomb lattice get shifted 
away and achieve different potential values. A single 
band tight-binding description will be appropriate for 
Vg when the wells are deep enough and sufficiently well 
separated in energy. The recoil energy associated to the 
"magic" honeycomb potential is here Er = EmCos^d, 
where E„i = h'^k'^J{2M) is the "magic" recoil energy, 
M being the mass of Yb atoms. Under the current ex- 
perimental configuration, one has E^ = AEam,/'^ where 
Earn = ^^fcam/(2-M) is the "anti-magic" recoil energy. 
Direct tunneling between theses two wells is then sup- 




FIG. 4. [Color online] Plot of Vg{x, y) as a function of x when 
y = l/\/3 for various Vi/Vq ratios. Black thin line: "magic" 
potential for V\ — 0. Blue thick line: Vi — Vo- Red dashed 
line: Vi — 3Vo. Green dotted line: Vi = 5Vb. 

pressed when Vi ^ Vo ^ Eam (we take Vi = 3Vo in 
our subsequent calculations). One can also show that the 
harmonic approximation around the global minima of Vg 
is no longer isotropic but features two different trapping 
frequencies along Ox and Oy. In turn, the correspond- 
ing Wannier functions Wg{r — Vg), centered at the \g)- 
sublattice sites rg, will reflect this spatial asymmetry. 
For Vi — 3Vb, the harmonic anisotropy is small and the 
corresponding harmonic lengths differ by 10% only. Ob- 
viously, the same conclusions apply to Vg and its asso- 
ciated Wannier functions We{r — re) centered at the |e)- 
sublattice sites r^- Furthermore, because of the inversion 
property between the two optical potentials one can infer 
the interesting property: 

w,i-r)^wg{v). (10) 

The superposition of these two independent regular 
triangular state- |e) and state- jg) sublattices defines a 
new honeycomb structure and Fig. [5] shows its primi- 
tive Wigner-Seitz cell. As one can see, contrary to the 
"magic" potential, the resulting new honeycomb struc- 
ture is no longer regular as the two sublattices are no 
longer shifted by a along Ox but by a' « 0.88a. The new 
nearest-neighbor vectors c'^ {a = 1, 2, 3) no longer add up 
to zero but, because the underlying Bravais lattice is still 
triangular, they still verify ai = c[ — and a2 = c!-^ — Cg . 
As a consequence, the overlap of the Wannier functions 
Wg and Wf, along the horizontal link will be different from 
their overlaps along the two other links (these being iden- 
tical), a fact which has its importance when coupling the 
two sublattices. 



III. UNIFORM FLUX CONFIGURATION 

So far we have been able to produce a honeycomb lat- 
tice where each of its sublattice traps atoms of a given 
internal state. However these sublattices are still uncou- 
pled as atoms cannot yet flip their internal state and hop. 
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FIG. 5. The Wigner-Seitz cell of the new honeycomb struc- 
ture created by the superposition of the two triangular state- 
dependent trapping potentials. This unit cell is no longer reg- 
ular and we have a' = 0.88a, a" = 1.07a. The proportions on 
the figure have been exaggerated. 



For this one needs to expose the atoms to Raman lasers 
which, by resonantly coupling the two internal states of 
the atoms, will induce these hops and thus couple the 
two sublattices. As will be explained below and in Sec- 
tion IIV( uniform or staggered synthetic magnetic fields 
can then be implemented for a suitable choice of the Ra- 
man lasers, and atoms traveling around a unit cell will 
pick up a nonzero phase. The net flux per cell can be 
made of the order of one quantum flux, thereby reaching 
the strong field regime. 



A. Raman-induced hopping 

The Raman coupling between the two honeycomb sub- 
lattices makes an atom with internal state \g) at site 
hop to a site at while flipping its internal state to |e) 
(and vice versa). For a plane- wave Raman laser field with 
wavevector q, the associated (complex) hopping ampli- 
tude reads: 



Jg_,e = ^ j fi2rw*(r-re)e"i ""u;g(r- 



where 17 is the Raman laser Rabi frequency. For the re- 
verse hopping process, one has simply Je->g = Jg^e- 
Most importantly, if one now considers the hopping am- 
plitude associated to Bravais-translated sites r'g = re -|- R 
and r' — To+R, where R is a Bravais lattice vector, then: 



J, 



g'- 



(12) 



Another interesting property can be found by defining 
Rge = (re + rg)/2 and r^e = re - r,,. Then: 



Ja 



(13) 



Se X iq.r 



Wo 



2 ' 



(14) 



Using now Eq. ([TU)l . it is easy to show that Ig^ = I*^ 
meaning that Ig^ is in fact real. Assuming the Raman 
laser Rabi frequency to be real, the phase of the hopping 
amplitude is then simply given by: 



(15) 



In the following, we will assume that the overlap 
between the Wannier functions is only significant for 
nearest-neighbor sites. This will be the case when Vi > 
Vo >• Earn- As Stated earlier, this overlap, and thus Ige, 
will be link-dependent as the new honeycomb lattice is no 
longer regular. We will thus neglect direct tunneling or 
second-order Raman-induced tunneling within each sub- 
lattice and only consider nearest-neighbor hopping be- 
tween the two sublattices. 
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FIG. 6. [Color online] Top figure; sketch of the hexagonal pla- 
quette alternating sites where sta.te-\g} atoms are trapped (A, 
C and E, open circles) and state-|e) atoms are trapped (B, D 
and F, red-filled circles). The total phase accumulated per cell 
is calculated for atoms hopping clockwise (arrows). Bottom 
figure: Energy diagram along the loop and the correspond- 
ing Raman transitions. All Raman transitions occur here at 
the angular frequency separation loq between the two hyper- 
fine state. As a consequence one single Raman laser beam is 
enough to address all sites. 
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B. One single Raman laser field is not enough 

Unfortunately, the simplest scheme where one uses a 
single Raman laser beam does not induce any global 
phase around a Wigner-Seitz plaquette. Indeed, consider- 
ing the situation depicted in Fig. [BJ one can identify pairs 
of hopping amplitudes related by a Bravais translation: 



T — o*qa2 T* 



(16a) 
(16b) 
(16c) 



It is then easily seen that the total phase accumulated 
around the Wigner-Seitz plaquette trivially cancels out 
and thus this Raman scheme alone fails to produce an 
artificial magnetic field. 



C. More is different 

The reason why the single Raman beam scheme fails 
is because all links are on an equal footing. To cure this 
problem, we need to consider a slightly more involved Ra- 
man laser scheme. To this end we introduce two counter- 
propagating laser beams with wavelength A producing a 
standing-wave pattern W{y) having the periodicity of Vg 
and Ve along Oy. This imposes A — 2-\/3a — 2Aam/"\/3 ~ 
1.29/im. These laser fields generate different optical po- 
tentials for state \g) and |e) (in units of Aa.,„/2): 
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(17) 



FIG. 7. [Color online] The left part of the figure shows the 
contour plot of the trapping potential Vg{Y). The right part 
of the figure shows the weak perturbing standing-wave po- 
tential which is added along Oy to lift the energy degener- 
acy between the potential wells. The standing-wave potential 
is maximum for y = and its period exactly matches \/3a 
where a is the "magic" honeycomb lattice constant. Every 
well in every other horizontal row is lifted up. Dark blue let- 
ters "G " mark the wells which are lifted up whereas blue 
letters "g" mark the unaffected ones. Accordingly, the dark 
red letters "E" and red letters "e" mark the respective po- 
sitions of the lifted-up and unaffected wells for the trapping 
potential Ve(r). Both sta.te-\g) and state-|e) lattices are rect- 
angular with a two-point unit cell (blue dashed-dotted and 
red dashed rectangles respectively), the corresponding Bra- 
vais vectors being (ai-|-a2) and (a2 — ai). The resulting global 
honeycomb structure exhibits two different vertical alternat- 
ing strips made either of the clockwise "gEGeGE" cell (dotted 
hexagon) or the "GegEge" one (solid hexagon). 



where we have asumed that the phases of these additional 
laser beams are fixed in such a way that the maxima of 
Wg^e along Oy coincide with those of Vm when x — Q. 
The new total potentials read: 



(18) 



Because of the chosen periodicity along Oy and choice of 
phase, the potential value is lifted in every other horizon- 
tal row of sites, the potential energy increase being SAg 
or h/S.e depending on the trapped internal state, see Fig. 

m 

In the following we will assume that hAg^e ^ Vb so 
that the net (perturbative) effect of Wg,e(r) is simply 
to lift the potential energy without modifying the origi- 
nal Wannier functions ■u;e,g(r)- However, because of these 
energy shifts, the honeycomb structure now features two 
different cells and exhibits vertical "strips" of these cells 
alternating along Ox, see Fig. [71 One now needs four 
Raman beams to address all lattice links and activate 
hopping between all neighboring sites, see Fig. [5] for one 
Wigner-Seitz plaquette and Fig.|9]for its neighboring one. 
Their respective angular frequencies are a;+ = loq + Ag, 
oj^ = ujQ — Ag, u) — LdQ + Ae — Ag and ujq and their re- 
spective wavevectors are q+. q_, q and qo. For sake of 



simplicity, we will assume that all these Raman beams 
propagate in the {Ox,Oy) plane. Since the light-shifts 
are small compared to the transition angular frequency, 
Ag e <C wqj if is legitimate to neglect the tiny changes in 
wavevector length so that we can consider that all these 
wavevectors have the same norm ko — 2tt/Xq. 

Then, for the Wigner-Seitz cell of Fig. [51 the cumulated 
total clockwise phase is simply: 



$ = (as -ai) 



while the one for that of Fig. [HI is: 



2q 



(a2 - ai) 



q+ 



2qo 



(19) 



(20) 



Of course we recover the fact that the phases cancel as 
they should when all wave vectors are identical. The in- 
duced artificial magnetic fiux will be uniform provided 
$ = <&', i. e. provided (q — qo) • {a.2 — ai) = 0. This will 
be the case if qo and q are most simply chosen parallel. 
If we now further choose q_(. and q_ to be parallel and 
anti-parallel to qo and q, the total cumulated phase then 
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reads <I> = 2-Ka with: 



a — — ;= — COS On. 



(21) 



where 9u is the angle of q_|_ with respect to axis Oy. The 
maximum value for a is 2Xam/{V3Xo) ~ 2.24, mean- 
ing that this scheme can provide a bit more than two 
flux quantum per cell. For comparison the correspond- 
ing magnetic field giving rise to the same maximum flux 
would he B = SttH/ {eXoXam) ~ 256G. To obtain the 
same flux in graphene samples, as the graphene lattice 
constant is very small (of the order of O.lSnm), one would 
need a magnetic field about 6 x 10® larger. 




FIG. 8. [Color online] Top figure: clockwise "gEGeGE" pla- 
quette. Open and grey-filled circles mark the unafi^ected and 
lifted-up |gi)-sites whereas red-filled and dark-red-filled circles 
mark the unaffected and lifted-up je)-sites. Bottom figure: En- 
ergy diagram for the Raman transitions around the plaquette. 
Because the wells are lifted differently, three different angular 
frequencies, u, 10+ and uj-, are now needed to activate the 
hops around this plaquette. 



IV. STAGGERED FLUX CONFIGURATION 
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FIG. 9. [Color online] Top figure: clockwise "GegEge" pla- 
quette. Open and grey-filled circles mark the unaffected and 
lifted-up |gi)-sites whereas red-filled and dark-red-filled circles 
mark the unaffected and lifted-up |e)-sites. Bottom figure: En- 
ergy diagram for the Raman transitions around the plaquette. 
Because the wells are lifted differently, three different angular 
frequencies, ujq, iJ+ and uj-, are now needed to activate the 
hops around this plaquette. 



by simply taking the q+ and q_ wave vectors to be 
along axis Ox and the qg and q wave vectors to be anti- 
parallel, the accumulated phases in the two different type 
of Wigner-Seitz cell would now be opposite: 



$' = -$ = (a2 - ai) • q = na, 



(22) 



where 6*^ in Eq. (|2f|) is now the angle of q with respect 
to Oy. We get alternating vertical stripes with opposite 
fluxes, the maximum value reached being half the uni- 
form one, i. e. about l.f 2 flux quanta. Hence, the config- 
uration proposed in Section IIIII proves quite versatile as 
it allows to easily switch from the uniform to the stag- 
gered magnetic flux cases by only changing the direction 
of propagation of the Raman beams. 



A. Adapting the previous configuration 



B. A novel configuration 



An artiflcial staggered magnetic flux can be easily cre- 
ated using the results of the previous Section. Indeed, 



However, for sake of completeness, we would like to 
propose here another scheme, similar to the one proposed 
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FIG. 10. [Color online] Staggered magnetic field case. Con- 
tour plot of the "magic" honeycomb lattice 'K(i(r) and the 
"anti-magic" standing- wave potential VamM with spatial pe- 
riod 3a, a being the "magic" honeycomb lattice constant. The 
blue letters "g" mark the positions of the state-lg) sites for the 
trapping potential Vg = + Vam while red letters "e" mark 
the positions of the state- |e) sites for the trapping potential 
= — Vam- Both sta.te-\g) and state- |e) lattices are rect- 
angular but with a two-point unit cell (blue dashed-dotted 
and red dashed rectangles respectively), the corresponding 
Bravais vectors being (ai +3.2) and (a2 — ai). The Raman 
laser beams only connect horizontal nearest-neighbor state- 
\g) and state- |e) sites. The resulting honeycomb structure ex- 
hibits two different alternating vertical strips made either of 
the clockwise cell "eeggge" (dotted hexagon) or the "ggeeg" 
one (solid hexagon). 



in [1^ , where the final lattice configuration alternate zig- 
zag vertical rows of |g)-state minima with zig-zag vertical 
rows of |e)-state minima. To achieve this, we simply start 
with a "magic" honeycomb potential where all beams 
are now coplanar with respective consecutive angles of 
120° like done in js^l- With respect to the previous sec- 
tion, this configuration is obtained by setting 6 — 0^ see 
Eq. ([T]). Using now Am/2 as the space unit in this Section, 
we get V^^ir) = V^,'[3 + 2/;„(r)] with: 
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FIG. 11. [Color online] Staggered magnetic field case. Plot 
of the trapping potential Vg{r) as a function of x when 
y — y/3a/2 for various V{/Vq ratios. Black thin line: "magic" 
honeycomb potential (Vi' = 0). Blue thick line: V{ — Vq. Red 
dashed line: Vl = SFq- Green dotted line: V/ = SVq. 
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FIG. 12. [Color online] Staggered magnetic field case. Top fig- 
ure: clockwise "eeggge" plaquette. The red-filled circles mark 
the state-je) sites while open circles mark the state-jg) sites. 
Bottom figure: Energy diagram for the Raman transitions 
around the plaquette. A single Raman laser beam at frequency 
(jjQ is required here to activate hopping along the horizontal 
"e-g" links. Hopping along the vertical direction is due to di- 
rect tunneling through a potential barrier. 



This "magic" honeycomb lattice constant is a = 4/(3\/3) 
in dimensionless units (a « 0.29/im in full units). 

The "anti- magic" potential is now produced with two 
laser beams counter-propagating along Ox, their off- 
plane elevation angle being 9'. In dimensionless units. 
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the "anti- magic" potential reads 

K-„,(r)^F/ ^ + ^°^^^;-^^"^°» , (24) 

where b' = b'e^ — Att cos 6' / Xam^x- We now request the 
period of the "anti-magic" standing- wave to match 35, and 
the horizontal shift to be Xq — 3a/ 4. This imposes the 
off-plane angle to satisfy 



cosW = 



4A„ 



giving 9' « 50.35°. One then finds: 



V' (r) = V.' 

am \'- / ^ 1 



,1 



i{V3ttx/2) 



(25) 



(26) 



and V'g,e(r) = V'^(r) ± K'^(r). With this set of parame- 
ters, one finds that b' = (bi + b2)/2. In other words the 
Bravais lattice associated to the potential Vg and Ve is 
no longer triangular. It turns out to be rectangular with 
reciprocal lattice vectors b' along Ox and b" = \/3b'ey 
along Oy respectively, the Wigner-Seitz cell of the lat- 
tice of global minima being a two-point cell, see Fig. (TU] 
As Vg and Ve are still related by inversion and reflection 
about Oy, the inversion symmetry relating their associ- 
ated Wannier functions remains valid. Figure [TT] shows 
a plot of Vg when y = \/3a/2 for various ratios V{/Vq. 
A single band description for each rectangular lattice is 
appropriate when the condition V{ > Vq ^ E„i is met, 
Em being the "magic" recoil energy. Typically a ratio 
V{/Vq = 3 or larger is required. 

To couple the two shifted independent rectangular \g)- 
state and |e)-state sublattices, it is then sufhcient to use 
a single in-plane Raman laser at frequency wq and wave 
vector ko. We assume that the Raman laser only cou- 
ples Yb atoms (with different internal states) which are 
trapped in nearest-neighbor sites. As Fig. [TU] shows, the 
Raman-induced hops only occurs along Ox, between ver- 
tical zig-zag rows. The hopping along the zig-zag rows 
relies on direct tunneling through the potential barriers. 
As a whole, the coupled system displays a honeycomb 
structure. Each vertical row is built by repeated tiling 
of the same plaquette. There are two different kind of 
plaquettes, obtained from each other by reflection about 
their middle vertical axis. The rows with different plaque- 
ttes are alternating along Ox. A schematic picture of the 
hops around the clockwise "eeggge" plaquette is shown 
in Fig. [121 The situation for the other "ggeeeg" plaquette 
is simply obtained from it by mirroring the sites through 
the middle vertical axis and changing the energy diagram 
accordingly. 

By an argument similar to the case of the uniform syn- 
thetic magnetic flux, the phase of the hopping amplitude 
between a |g)-state site and its neighboring |e)-state site 
is simply given by ko • Rge, where is the middle point of 
the (horizontal) connecting link, the phase for the reverse 
hopping process being the opposite. As the hopping am- 
plitude along the vertical zig-zag chains is real, the total 



phase accumulated around the cell shown in Fig. [TO] is 
simply 5- = 27ra' = (47rAm/3Ao) cos0^, where is the 
angle of the Raman wave vector with axis Oy. The phase 
for the other plaquette is — ^E". The maximum value ob- 
tained for a' with this scheme is (2A„i/3Ao) ~ 0.876, a bit 
less than a flux quantum. Here again we get alternating 
vertical stripes with opposite fluxes. 



V. HARPER MODEL 

In this section we would like to give the single-band 
Harper model [2l| . [2^ . [4^ describing the previous config- 
urations and discuss some orders of magnitude and lim- 
itations. We will restrict ourselves to the configuration 
obtained in Section IIIII as it can as well describe a uni- 
form or a staggered synthetic magnetic flux applied to a 
honeycomb lattice. 

To find the Harper Hamiltonian associated to the con- 
sidered optical potential configuration, one has first to 
identify its unit cell C and the corresponding Bravais lat- 
tice B. Then the Harper Hamiltonian is simply written 



-^Harper — ^ ^ , 



(27) 



where i7„ is the hop operator acting on the cell 0„ ob- 
tained by translation of C along B. In turn, one can write: 

Hn = Tn + Tl (28) 

where T„ describes hops among sites within 0„ but along 
prescribed directions. 



A. Simplest uniform flux configuration 

This is the configuration obtained for q = qo = ^ q+ = 
— q_. In this case, it is easy to see that the lattice is 
obtained by repeated tiling of the unit cell displayed in 
Fig. [TOl the relevant Bravais lattice B being triangular 
and spanned by ai and a2. We take the origin of the 
cell at the state- jg) site located at Tn. We denote by 
the corresponding annihilation operator and by Can the 
annihilation operators of a state- |e) atom located at sites 
Tan = Tn — (a = 1, 2, 3). The hop operator is: 



(29) 



where Jq is the (real) hopping strength along the link 
vector -c^ and (j)an = £a q- (r„ -</2) (^2 = £3 = -£i = 
1). As obtained from Section iHll the total anti-clockwise 
phase per unit cell is: 



$5 = 2q- (a2 - ai). 



(30) 
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At this point one may wonder which gauge potential A 
and which magnetic field B = V x A would give rise 
to the same Harper model for electrons in the graphene 
lattice. Using the gauge potential Al — —Byex (Landau 
gauge), we write A = A^ + Aq. The vector potential 
generates a uniform magnetic field perpendicular to the 
honeycomb lattice and we choose its strength B such that 
it gives rise to the same flux per unit cell as A. It is then 
easy to show that the vector potential Aq, which does 
give rise to a nonzero magnetic field, nevertheless gives 
rise to zero flux around any closed loop. In the Harper 
model the particles move along the lattice links, it means 
that the contribution of Aq can in fact be locally gauged 
away. Indeed, picking up some lattice site ro as the origin, 
the phase = /^J* Aq • dr is in fact path independent 
(as long as the path is taken on the honeycomb lattice) 
and only depend on the end point r„. This proves that 
the Harper model given by Eqs d^T]), (QSl) and (gl]) can 
be in fact obtained with a uniform magnetic field. 
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FIG. 13. [Color online] The cell used to derive the Harper 
Hamiltonian for the simplest uniform magnetic flux configu- 
ration. Red circles: state-|e) sites. Open circle: sta.te-\g) site. 
The arrows on the links indicate the hop directions consid- 
ered to write down the hop operator T„, see Eq. (|28p . The 
Raman wave vectors implied in the hops are indicated along 
each link. The hop operator operates on the same cell 
but with reversed hop directions and thus complex conjugate 
hopping amplitudes. The Harper Hamiltonian is obtained by 
superposing the hop operators obtained by repeated transla- 
tions along the Bravais lattice B spanned by vectors ai and 
a2, see Eq. (f27)l . 

We would like to give now a somewhat simpler ex- 
pression for the hop operator. Performing the local C/(l) 
gauge transformation: 

gn exp(-iq • r„/2) 

eln exp(-iq- ra„/2)eT„, 

one gets: 

T„ = Jie"^ e|„5„ -I- J2 ej„5„ -I- J3 el^gn (32) 

where (/? = q • (r„ + ri„) is the Raman phase evaluated 
at the mid-point along the link vector —c'l. One can eas- 
ily check that the total anti-clockwise phase picked up 



around a hexagonal plaquette is again $3. This local 
U{1) gauge transformation amounts to choose locally the 
phase of the Wannier functions. 

Because the link vectors and C3 have same length, 
the tunneling amplitudes J2 and J3 are the same as long 
as the corresponding Raman Rabi frequencies are equal. 
However, as c'l has a different magnitude, in principle Ji 
is different from the two others. The mismatch could be 
compensated for by a fine-tuning of the corresponding 
Raman Rabi frequency though this might not be easy in 
practice. In another context, it is known that hopping 
imbalance can have a strong impact on the physical phe- 
nomena under study (43| . 

Forgetting about this hopping imbalance, one can give 
a rough order of magnitude of the hopping amplitude 
J = {hn/2)Ieg, see Eqs. and ([T?l) . Using a harmonic 
approximation, and further neglecting the Wannier func- 
tions anisotropy which is anyway small, a simple calcu- 
lation performed at Vi — 3Vb shows that: 

nn 

J~ — exp[-(1.587/;ie + 0.323;ie)], (33) 

where the effective Planck's constant is — y/'2Eam/VQ 
and Earn is the anti-magic recoil energy [Eam/h — 
900Hz). For hD. — 2Eam and Vq = IGEam, one gets 



B. General case 

In the general case where all Raman wave vectors are 
different, the unit cell must now be enlarged and it is 
given in Fig. 1141 Writing down the hop operator T„ is 
straightforward but rather tedious as now many links are 
required. We leave it as an exercise for the reader. The 
noticeable point is that the relevant Bravais lattice B is 
now rectangular and spanned by the vectors ai + a2 and 
a2 - ai . 

The above Harper model can be readily extended when 
interactions come into play. This could be experimen- 
tally realized by loading the lattice with bosonic atoms 
in the presence of on-site repulsive interactions. In this 
case one could target the fractional quantum Hall ef- 
fect and the Laughlin state at filling fraction i' = 1/3, 
and more generally highly-correlated quantum liquids. 
By loading the lattice with fermionic atoms, one could 
target the relativistic quantum Hall effect as evidenced 
in graphene samples at low fluxes per cell. It is worth 
mentioning that contrary to ^^*Yb which has zero to- 
tal spin, the fermionic Yb isotopes have a nuclear spin. 
They are thus multi-level systems and one can then think 
of designing more elaborate configurations to mimic non- 
Abelian gauge fields. One could also imagine loading the 
lattice with Bose- Fermi mixtures like ^'''^Yb— ^"^"^Yb jill 
or like ^''^Yb— ^^"'Yb where next- nearest- neighbor inter- 
actions of the order of six times the tunneling energy 
have been reported at zero magnetic field |45|. The case 
of Fermi-Fermi mixtures raise tantalazing questions on 
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FIG. 14. [Color online] The cell used to to derive the Harper 
Hamiltonian in the general case. Red circles: state-|e) sites. 
Dark red circles: lifted-up state-|e) sites. Open circles: state- 
ly) sites. Grey circles: lifted-up state-\g) sites. The arrows on 
the links indicate the hop directions considered to write down 
the hop operator r„, see Eq. (|28p . The Raman wave vectors 
implied in the hops are indicated along each link. The hop 
operator operates on the same cell but with reversed hop 
directions and thus complex conjugate hopping amplitudes. 
The Harper Hamiltonian is obtained by superposing the hop 
operators obtained by repeated translations along the rectan- 
gular Bravais lattice B spanned by vectors ai + a2 and a2 — ai , 
see Eq. P7)) . 

spinor superfluidity. For instance, the ^^^Yb— ^^^Yb mix- 
ture realizes a system with SU{2) x SU{Q) symmetry 



VI. CONCLUSION 

In this paper we have proposed experimental set-ups 
realizing Abelian gauge fields acting on Yb atoms moving 



in a honeycomb lattice and giving rise to uniform or stag- 
gered synthetic magnetic fluxes. A net flux per unit cell 
of one quantum flux can be easily reached, thereby real- 
izing the cold atom analogue of the Harper model. Dif- 
ferent phenomena could be experimentally studied with 
these configurations, ranging from the relativistic to the 
fractional quantum Hall effects. A possible extension of 
this work would be to study the role of the honeycomb 
lattice distortion which lead to hopping strength imbal- 
ance. As far as we know the impact of this imbalance on 
the Harper model is largely unexplored but would be of 
great experimental relevance. 
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